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The Belinkskii, Khalatnikov and Lifshitz conjecture says that as one approaches 
space-Hke singularities in general relativity, 'time derivatives dominate over spatial 
0^ , derivatives' so that the dynamics at any spatial point is well captured by an ordinary 

\ differential equation. By now considerable evidence has accumulated in favor of these 

CN \ ideas. Starting with a Hamiltonian framework, we provide a formulation of this 

^ ■ conjecture in terms of variables that are tailored to non-perturbative quantization. 

^ Our formulation serves as a first step in the analysis of the fate of generic space-like 

\^ , singularities in loop quantum gravity. 
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I. INTRODUCTION 

>■ I In 1970, Belinskii, Khalatnikov and Lifshitz (BKL) made a conjecture which, if true, 

would shed considerable light on the nature of space-like singularities in general relativity [H . 
They suggested that as one approaches these singularities, time derivatives would dominate 
over spatial derivatives, implying that the asymptotic dynamics would be well described by 
an ordinary differential equation. At first this claim seems astonishingly strong. However, 
over the years considerable analytic and numerical evidence has accumulated in its favor 
O ■ see, e.g., 0, i, 3). Together, these results suggest the following scenario: i) Geometry at 

any spatial point is well described by the Bianchi I metric for long stretches of time; ii) 
From time to time the parameters pi characterizing the Bianchi I metric undergo a (Bianchi 
II) transition, rapidly settling to new values and, iii) The spatial gradients can grow but 
they do so on 'small sets' where spikes develop.^ In addition, except for a scalar field or a 
stiff fluid, sources do not play much of a role in dynamics; 'matter doesn't matter' close to 
the singularity. 

It is tempting to use this scenario as a starting point in the analysis of what happens to 
these classical singularities in the quantum theory. In particular, Garfinkle [5] has suggested 
that understanding of the quantum behavior of the Bianchi I model would shed considerable 
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^ As is common in the recent literature, we use the convention that the singularity occurs at f = — oo. So 
in general there are infinitely many transitions till one reaches the singularity. However, general relativity 
cannot be trusted once the curvature acquires Planck scale and until this epoch there are only a finite 
number of transitions and spikes at any spatial point. 
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light on the fate of space-hke singularities in quantum gravity. It is now known in loop 
quantum cosmology that the Bianchi I singularity is naturally resolved because of quantum 
geometry effects [1, 0]. This suggests that there may well be a general result which says 
that all space-like singularities of the classical theory are naturally resolved in loop quantum 
gravity. 

However, it is difficult to test this idea using the current formulations [sl, l9|] of the BKL 
conjecture. For, these formulations are motivated by the theory of partial differential equa- 
tions rather than by Hamiltonian or quantum considerations. The basic variables, for ex- 
ample, are chosen to simplify functional analysis and numerics. While this is a natural 
strategy within classical general relativity, these variables are not well suited for quantum 
theory. For example, one often uses the variables T,ab, Nab, introduced by Uggla, van Elst, 
Wainwright and Ellis (UEWE), which are obtained by dividing certain geometric fields by 
the trace K of the extrinsic curvature P]. It is difficult to write down the corresponding 
operators in the quantum theory particularly because of the factors. So, the question 
arises: Is there a formulation of the BKL conjecture in terms of variables that are well suited 
for non-perturbative quantization? The answer turns out to be in the affirmative. The goal 
of this communication is to present a succinct summary of this formulation. As we will see, 
the framework has some features which make it appealing also in the classical theory. 

We will focus on vacuum general relativity because 'matter does not matter'. Details 
for the scalar field — which does matter — have been worked out and key differences will be 
summarized in section IVl 



II. CONVENIENT VARIABLES 

We will consider space-times (*^M, '^Qab) with = M x M, where M is a compact 3- 
manifold without boundary. Let us choose as our gravitational variables pairs {Ef, KD 
of fields on M, where a,b,c. . . are tensor indices and i,j,k. . . are S0(3) internal indices 
(which can be freely raised and lowered using the Cartan Killing metric on so (3)). E°- 
represents an orthonormal triad on M (with density weight 1) and -K"* represents extrinsic 
curvature 'on shell'. These are canonically conjugate on the gravitational phase space [ll|, 
Q: {E-{x), Kliy)} =5^5i5%x,y). 

The density weighted triad E°- determines a positive definite metric qab on M via EfE^'^ = 
qq°'^ where q is the determinant of qab- The standard extrinsic curvature Kab is given by 
E'"Kab = y/qKl^. The scalar, vector and Gauss constraints of vacuum general relativity are 
given by 

S:=-qR-2E[,E]^KlKl^Q, Va := ^ D^a [KI^E]) ^ Q , G" := u^'EJKJ ^ , (2.1) 

where D and R denote the derivative operator and the Ricci scalar of qab- As usual, the 
Hamiltonian generating dynamics is just a linear combination of these constraints. The 
triad Ef determines a unique S0(3) connection FJ^ through DaE^ + eij'^TlE^ = 0. One can 
define another SO (3) connection via := F^ — 7 i^* where 7 is called Barbero-Immirzi 
parameter. Loop quantum gravity is based on the canonical pair (A^ E^). This pair provides 
the point of departure for our formulation of the BKL conjecture. 

In the classical analysis we now wish to undertake, it is simpler to use all three fields, 
F^, i^*, E"-, although F^ is determined by E"-. The key idea can be then summarized as 
follows. The accumulated results to date suggest that the metric qab becomes degenerate 
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at the singularity, whence its determinant q vanishes there. (For example in the Bianchi I 
model, in terms of the commonly used proper time t, the metric is given by ds^ = —dt^ + 
Ylii't'^^^dxl and the singularity occurs at t = 0. Since = 1, we have q = t^.) Therefore 
one might expect that fields which are rescaled by appropriate powers of q would remain well 
behaved at the singularity. Similarly, while the covariant spatial derivatives Daf of fields 
/ may not be sub-dominant compared to time-derivatives, derivatives Dif := E^Daf are 
more likely to be, because of the ^ factor in the density weighted triad Ef. This strategy 
is similar to that used in a more common formulation of the conjecture j2|, Q, [sl where, as 
mentioned in section [H one divides geometric fields by the trace K of the extrinsic curvature 
which is expected to diverge at the singularity. The relation between the two strategies is 
discussed in section IIIII 

These motivations lead us to introduce the following basic variables: 

^=EtVi-Eirl5l and ^= E^Ki - E^KI 51 . (2.2) 

It turns out that constraints can be re-expressed entirely in terms of Cij,Pij and their 
derivatives: 

S := 2e'''D,{C,k) + 4%]C[^^] + QjC^' - + P,,P^' " ~ (2-3) 

:= -2D,P^ - 2e,klP^' C' + t:jk{2P^'Ci^ - CP''') ^ (2.4) 
Qk ._ ^ijkp..^Q_ (2.5) 

Next, let us consider evolution equations, i.e., the Hamiltonian flow generated by the con- 
straints. For simplicity, in this brief communication we will set the shift to zero. As in loop 
quantum gravity, our lapse N will be a scalar density of weight -1. Then the time evolution 
of our basic triplet {Cij, Pij, Di) is governed by: 

= \ e/'D,{N{2Pu - 6uP)) + N [2C(,|fc|P^-) + 2C[kj]P\ - PC,,] (2.6) 

Pi, = -ej^'DkiNCu) + ^e,,'Dk{NC) - e''"'DUNCki)S,, (2.7) 

+2e,fc™C[,fc] DmN + {D,Dj - dijO^Dk) N 
+N[-2C(^ik)C''j + CCij - 2C^''^^C[ki] Sij] 

ft 

D^Sn = - [AiVP]s„ - NP/D,Sn (2.8) 

where, in the last equation, Sn is any scalar density of weight n. These equations can be used 
as follows. On an initial slice, we construct {Cij, Pij, Di) from a pair {E'^,KD of canonical 
variables. But then we can deal exclusively with the triplet {Cij, Pij, Di). The pair {E!l, KD 
satisfies constraints if and only if the triplet satisfies fl2.3l) - fl2.5p . Given such a triplet, we 
can evolve it using fl2.6p - fl2.8l) . without having to refer hack to the original canonical pair 

iE:,K). 

These two sets of equations have some interesting unforeseen features. First, the basic 
triplet {Cij, Pij, Di) has only internal indices: our basic fields are scalars (with density weight 
1). It would be of considerable interest to investigate if this fact provides new insights into 



the dynamics of 3+1 dimensional gravity [l3j. Second, these equations do not refer to the 
triad Ef. Suppose we begin at an initial time where Cij is derived from an E^. Then these 
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constraint and evolution equations ensure that Cij is derivable from a triad at all times. 
Furthermore, we can easily construct that triad directly from a solution {Cij,Pij) to these 
equations: first solve fl2.6p - fl2.8p and then simply integrate the ordinary differential equation 

= -NP/E^ (2.9) 

at the end. Third, the structure of the constraint and evolution equations in terms of 
{Cij, Pij, Di) is remarkably simple since only low order polynomials of these variables are 
involved. Finally, thanks to our rescaling by ^/q, our basic triplet Cij, Pij, Di (as well as Ef) 
is expected to have a well behaved limit at the singularity. (This expectation is borne out 
in the Bianchi I and II models.) Note also that our equations are meaningful even when 
the triad becomes degenerate. So, strictly (as in loop quantum gravity H, El]) we have 



a generalization of Einstein's equations. To summarize, we have found variables which are 
likely to remain finite at the singularity and rewritten Einstein's equations as a closed system 
of differential equations in terms of them. Therefore, this formulation may be particularly 
useful for proving global existence and uniqueness results . 

However, in the formulation given above, all variables carry density weights. Also, al- 
though the operator Di satisfies linearity and the Leibnitz rule, it is not a standard derivative 
operator: it has torsion and, more importantly, it changes the density weight by 1. These 
unconventional features may make these equations awkward to handle particularly for nu- 
merical work. In the detailed paper [13] we will provide an equivalent formulation involving 
only functions (with zero density weight) and a more familiar version of the operator Di. 



III. THE CONJECTURE 

The key step in any formulation of the BKL conjecture is to specify the basic variables, 
what one means by their 'spatial derivatives' (which are to be sub- dominant), and 'time 
derivatives' (which are to dominate). We will consider any smooth foliation Mt of an appro- 
priate portion of such that the space-like singularity of interest constitutes a (limiting) 
leaf. The time function t labeling our spatial slices is intertwined with the choice of lapse. 
We will assume that the density weighted lapse N admits a smooth limit as one approaches 
the singularity. Now, we are led by the intuition that the spatial metric qabif) becomes 
degenerate at the singularity. This implies that the lapse function N (with density weight 
zero), given by N := ^/qN , goes to zero, i.e., that the singularity is at t = — oo. 

In this set up, our basic variables will be {Cij, Pij) and the lapse N . By time derivatives, 
we will mean their Lie derivatives along the vector field t°- := Nn°' where is the unit normal 
to the foliation Mt. By spatial derivatives we will mean their Di derivatives. Since Di := 
E^Da, the notion does not depend on coordinates. Rather, it is tied directly to the physical 
triads and the covariant derivatives compatible with the metric. Then, the idea behind the 
conjecture is that, as one approaches the singularity, the spatial derivatives DiCjk, DiPjk 
of the basic fields should become negligible compared to the basic fields themselves (in 
particular) because of the y/g multiplier in the definition of E^. An immediate consequence 
is that the antisymmetric part of Cij is negligible [lo|, a fact that we will repeatedly use 
below. 

Thus, our formulation of the BKL conjecture is that, as one approaches the singularity, 
solutions to the Einstein's equations (2.3)-(2.8) are well approximated by solutions to the 
truncated system of equations obtained by ignoring the spatial derivatives. The truncated 
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constraints are purely algebraic equations: 

:= C,,C^' - + P,,P^' - ~ (3-1) 

l^f ^ := e,,fe(2P^'G'= - CP''') ^ (3.2) 
Gfr) := e*^'=P,,^0, (3.3) 

while the truncated evolution equations are ordinary differential equations: 

Q,- = iV [2Cfe(iP',)) - PC,,] and P^, = iV[-2CifcC^- + CC^,] . (3.4) 

Note that operation of truncation and hence the final truncated system depends crucially on 
one's choice of basic variables and notions of space and time derivatives. (For example, if we 
had used triads rather than Cij as basic variables, we would have been led to set Cij to zero 
in the truncation procedure, and truncation would have led us just to Bianchi I equations. 
The resulting BKL conjecture would have been manifestly falsej^ 

How does this formulation compare with that of UEWE 0, [1]? In that framework, one 
divides the geometrical fields by the trace of the extrinsic curvature K (which is expected to 
diverge at the singularity) while here we multiply them by the volume element ^Jq (which 
is expected to go to zero). There, the (scalar) lapse iV is such that NK admits a limit N. 
Therefore iV goes to zero and, as in our case, the singularity lies at t = — oo. The key scale 
invariant variables (A^jj,Ejj) which are expected to be well behaved at the singularity are 
related to our (Cjj, Ki^) via 

f^.. = Qp-^ and = -6P-^P(ij) + 25ij , or, (3.5) 

Cij = -^V^Nij and P(ij) = (S^j - 2%) . (3.6) 

Similarly the two sets of lapse fields are related simply by: N = NKy/q. Thus, although 
the motivations and the starting points of the two frameworks are quite different, the basic 
variables are closely related. From the viewpoint of differential equations, the two reduced 
systems would in essence be equivalent if K^/q admits a finite, nowhere vanishing limit at 
the singularity. This condition holds for Bianchi I models (and also Bianchi II which describe 
the transitions between Bianchi I epochs). An advantage of the {Cij, Pij) framework is that 
it is better adapted for non-perturbative quantization since it comes from the Hamiltonian 
framework underlying loop quantum gravity. 
Remarks: 

i) In our formulation we have allowed a large class of foliations. However, a closer examina- 
tion from the standpoint of differential equations may well lead to further restrictions. The 
inverse mean curvature foliations commonly used in conjunction with the UEWE framework 
appear to be well-suited also for our framework. 

ii) As a rule of thumb one often says that scale invariant scalars (such as Nij and Ej^) should 
have well-defined limits at the singularity. Those considerations can be extended to density 
weighted quantities used in this paper. The rule of thumb then is that a density of weight n 
has well-defined limit if it has scaling dimension 2n. Since Cij and Pij have density weight 
1 and scaling dimension 2, they should have a well defined limit at the singularity. 

iii) Consider the subspace of the full phase space on which DiCjk = and DiPjk = and 
demand that the lapse N satisfy DiN = 0. Then, a non-trivial result is that the Hamiltonian 
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vector field of full general relativity is tangential to this sub-space; conditions DiCjk = 
and DiPjk = are preserved by the full evolution equations fl2.6l) - fl2.8p . 

IV. THE BKL TRUNCATED HAMILTONIAN SYSTEM 

Let us now explore the BKL truncated system we were led to in section UTTl by focusing 
on the pair {Cij,Pij). We already know that Cij is symmetric due to the BKL truncation. 
The Gauss constraint implies that Pij is symmetric. Therefore, the basic variables of the 
truncated theory are pairs of symmetric fields {Cij, Pij) on M. 

Upon truncation, the full symplectic structure of general relativity yields the following 
Poisson brackets between these fields: 

{Cijix), Cki{y)}T = 0, and {Pij{x), PH{y)}T = 0, 

{Pij{x), CM{y)}T = {CkijS,)i + Ci(j5,)k){y) 5\x,y) . (4.1) 

These are subject to the truncated scalar and vector constraints (13.11) and (13. 2p . This system 
of constraints is of first class. As before, for simplicity, let us set the shift to zero and obtain 
the evolution equations by taking Poisson brackets with the scalar constraint. The result is: 

CiJ=N[2Cki^Pj)^ ~PCi,i and = -N [2CuCj'' - Cdj] . (4.2) 

This is the Hamiltonian flow generated on the truncated phase space by the truncated 
constraints. It agrees with the evolution equations (13.41) of section [TTll obtained by first using 
the full constraints on the full phase space and then truncating the full evolution equations. 
Thus, the truncation procedure is in complete harmony with the Hamiltonian framework. 

To explore the structure of the truncated theory, it is convenient to solve and gauge fix the 
vector constraint. Since our fields Cij and P^j are now symmetric, the vector constraint (13. 2p 
implies that, regarded as matrices, they commute, whence we can simultaneously diagonalize 
them. Finally, we can gauge-fix this constraint by demanding that Cij and Pij be diagonal. 
Then we arrive at the following description of the truncated phase space. It is coordinatized 
by the diagonal elements C/,P/ of Cij and Pij (with (/ = 1,2,3)). Their Poisson brackets 
are given by: 

{C,(x), CM}t = 0, {Pi{x), Pj{y)}T = 0, {Ci{x), P,(i/)}t = -25uCj5{x,y) . (4.3) 
There is a single (scalar) constraint: 

S^T){x) := \{^Cj{x)Y -^C^x) + ^(S Pj{x)Y -^Phx) ^0. (4.4) 

Thus, as expected, the truncated Hamiltonian system is ultra-local; dynamics at each spatial 
point is insensitive to what is happening elsewhere. Equations of motion are given by: 

Cj{x) = -NCj{x) (SPj(x) - 2Pi{x)) , Pi{x) = NCiix) (SCj(x) - 2Ci{x)) . (4.5) 

Note that, although the constraint is symmetric under interchange of Ci and Pj, there is 
a basic asymmetry in the symplectic structure which descends to the equations of motion. 
Dynamics leaves the sector of the phase space with Cj > invariant. Let us focus on this 
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sector and set Qj = — \n Cj/2. Then, (Q/, Pi) are canonically conjugate. Rewriting the con- 
straint (14.41) in terms of Q/, one immediately sees that what we have is particle dynamics in 
exponential potentials. By making a rigid rotation in the truncated phase space, (14.41) can be 
brought to the famihar Misner form S'(t) = -Pq^ + P^ + P^ +e-(^/^)'3o V{Q±). However, 
Oo) Q± are components of a connection while the Misner variables refer to components of 
the 3- metric. 




FIG. 1: Past evolution of the three Ci at a fixed point x, given by Eq ()4.5p . As time decreases, 
C2 starts growing rapidly near t = — 1, reaches a peak, and then goes to zero rapidly. Then near 
t = —23, C3 follows the same pattern and finally Ci does the same near t = —230. Their profiles are 
well described Eq (j4.6p . This growth and decay repeat ceaselessly as the singularity is approached 
at t = —00. In this simulation, the initial data at t = was Pi = —0.4742,^2 = —3.3586,^3 = 
-1.3096; Ci = 0.0746, C2 = 0.0040, C3 = 0.0070. 



Note that if C/(x) vanish initially at a point x in M, they remain zero throughout the 
evolution and Pi{x) remain constant. Then the dynamics of fields at that point x is that 
of the Bianchi I model. ^ (The P/ are linear combinations of the parameters Pi normally 
used to characterize the Bianchi I metric.) Dynamics of generic initial data is much more 
complicated than what one might expect from the deceptively simple form of the evolution 
equations (14.51) . To get a handle on the Bianchi II transitions, one can linearize (14.51) about 
a Bianchi I solution. Since we have assumed that the universe contracts as t decreases, it 
follows that all three P/ of the Bianchi I solution are negative. One finds that as one evolves 
backward in time two of the the three modes (C/, P/) are stable and one, with largest |P/(a;)|, 
is unstable. Let us suppose that |Pi(a;)| is the largest. Then the exact evolution of the stable 
modes 2 and 3 is well tracked by the linearized equations which imply C2{x, t) ?a ~ C^{x, t) 
and P2{x, t) ^ P2{x, to), P3{x, t) ~ P3(x, to)- However, for the unstable mode we have to use 



^ Note that since M is allowed to have any compact topology, it may not admit a global Bianchi I solution. 
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FIG. 2: Past evolution of the three Pj at a fixed point x (in the same simulation as in Figured]). 
Their profiles are well described by Eq (j4.7|) . Only one of the three Pj changes at a time, causing 
the solution to make a transition from one Bianchi I metric to another. The jump in each Pj is 
correlated in time with that in the corresponding Cj. 



the full equations (14.51) . Setting C2 = C3 = one obtains the solution: 

Ci{x, t) = ±2(/^ sech (2/^iV {t - r)), and (4.6) 
Pi{x,t) = P2 + P3 + 2v/i^tanh(2v/^iV(t-r)) (4.7) 

for some constant r, where the sign in the first equation is plus (minus) if Ci{x, to) is positive 
(negative). Thus, only one pair, namely {Ci{x,t), Pi{x,t)), evolves non-trivially. The form 
of solutions (14. 6 p implies that, as we evolve backward in time, there is a quick transition 
in which |Ci(x)| first increases, reaching a maximum value at t = r, and is then driven to 
zero, while Pi{x) is mapped to P2{x,to) + P3{x,to) + 2^y P2{x, tQ)P^{x, to). Simple algebra 
shows that this transition from one Bianchi I solution to another is given precisely by the 
well-known u-map. This derivation of the u-map does not require one to approximate the 
exponential potential in (14.41) by a rigid, perfectly reflecting wall. Also, our parametrization 
of the Bianchi-I solutions by Pj (rather than by the standard Pi) makes it easier to follow 
the Bianchi II transitions both analytically and numerically because in any transition only 
one of the three Pj changes (while all three pi change). Figures [1] and and [2] show these 
transitions in a numerical solution to the ODEs (14.51) . Here, an iterative Picard algorithm 
was used with an adaptive time step, chosen to ensure conservation of the scalar constraint. 
Such transitions in the truncated model mimic the situation in the full theory surprisingly 
well. (Compare, for example, our figures [1] and [2] with figures 5 and 4 in the second paper 
of i). 

So far we focused on the dynamics at a single spatial point x. Let us now consider a 
2-dimensional surface in the physical space on which one of the Cj, say Ci, vanishes and 
C2, C3 are small. Fix a point x on this surface. Then, Ci is positive to one side of the 
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2-surface and negative to the other. At a point ?/ in a neighborhood of x, as t decreases, 
Ci{y) first increases rapidly if it is initially positive and decreases rapidly if it is initially 
negative (following the sech(t — r) profile of fl4.6p ). This produces steep gradients at x for 
some time, which appear as spikes. However, as t decreases further, |Ci| goes to zero rapidly 
on both sides of the 2-surface for t < t. So the spike in Ci becomes dilute and disappears 
(but it will reoccur and disappear again repeatedly on further backward evolution). Pi on 
the other hand remains constant at x but decreases on either side of the 2-surface. As t 
decreases, the gradient of Pi keeps increasing and the spike sharpens. These transitions and 
spikes have been observed in numerous simulations of the truncated system 15|]. However, 



because the spatial gradients become large at spikes, the truncated system becomes a poor 
approximation and one has to consider the full system (as, e.g., in [16]). 

The analytical arguments given above explain all the qualitative features of the truncated 
dynamics if the initial phase space point lies in a neighborhood of a Bianchi I fixed point. But 
we do not have an analytical understanding of what happens away from this neighborhood. 
Numerical simulations show that even if one starts far away from pairs {Cj, Pj) corresponding 
to Bianchi I solutions, dynamics drives the system quickly to the Bianchi I sub-space. It 
would be very useful to have an analytical derivation of this phenomenon for our system 
( 14. 5 p without having to make further approximations. 



V. DISCUSSION 



We began with the Hamiltonian formulation of general relativity underlying loop quantum 
gravity where the basic fields are spatial triads E°- with density weight 1, spin connections 
they determine, and extrinsic curvatures K^. Based on the examples that have been studied 
analytically and numerically, the general expectation is that the determinant q of the spatial 
metric qat would become degenerate and the trace K of the extrinsic curvature would diverge 
at space-like singularities. One can therefore hope to obtain fields which remain well-defined 
at the singularity either by multiplying natural geometric fields by suitable powers of q or 
dividing them by suitable powers of K. In a commonly used UEWE framework 0, @], one 
chooses to divide by K. The resulting fields satisfy differential equations with desirable 
properties. However, because of the presence of K~^, in quantum theory it is difficult to 
introduce operators corresponding to the new fields. 

We adopted the complementary strategy of multiplying geometrical fields by ^g. In the 
Hamiltonian formulation with which we began, the basic field E"" is already obtained by 
multiplying the orthonormal triad by ^/q. One would therefore expect it to vanish at the 
singularity and this expectation is borne out in examples. Consequently, Ef also provides a 
natural avenue to construct additional fields needed in the BKL conjecture. Our variables 
Ci^ and Pi^ were obtained (modulo trace terms) simply by contacting the spatial indices of 

and Ki by E"-. Furthermore, because E^ vanishes in the limit, the operator Di := EfDa 
provides a convenient tool to express the notion of 'spatial derivatives which are to become 
subdominant' near the singularity. The main expectation is that asymptotically DiCjk and 
DiPjk would become 'negligible' relative to Cjk and Pjk- Now, in exact general relativity, 
time derivatives of Cij and Pij can be expressed in terms of their Dt derivatives, purely 
algebraic (and at most quadratic) combinations of Cij and Pij, the lapse N and its Di 
derivatives (see fl2.3l) - fl2.8l) ). Therefore, if in the limit the Di derivatives of the basic fields 
become negligible compared to the fields themselves, we are naturally led to conclude that 
time derivatives would dominate the spatial derivatives. This chain of argument led to our 
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formulation of the BKL conjecture. 

This rather simple idea depends on the fact that the structure of Einstein's equations 
has an interesting feature: as saw in section [Tll once the triplet Cij,Pij,Di is constructed 
from the triad Ef and the extrinsic curvature on an initial slice, the constraint and 
evolution equations can be expressed entirely in terms of the triplet. Given a solution to 
these equations, the spatial triad Ef (and hence the metric qat) can be recovered at the end 
simply by solving a total differential equation, (12. 9p . This is a surprising and potentially 
deep property of Einstein's equation. It provided key motivation for our formulation of 
the BKL conjecture and could well capture the essential reason behind the BKL behavior 
observed in examples and numerical simulations. Therefore, an appropriate quantization 
of the truncated system, e.g., a la loop quantum cosmology, could go a long way toward 
understanding the fate of generic space-like singularities in quantum gravity. 

Since the framework is developed systematically from a Hamiltonian theory, the BKL 
truncation naturally led to a truncated phase space. The specific truncation used has an 
important property: The truncated constraint and evolution equations on the truncated 
phase space coincide with the truncation of full equations on the full phase space. On the 
truncated phase space we could solve and gauge-fix the Gauss and vector constraints to 
obtain a simple Hamiltonian system. Solutions to this system have been explored both 
analytically [lo| and numerically [iH]- They exhibit the Bianchi I behavior, the Bianchi II 
transitions and spikes as in the analysis of symmetry reduced models and numerical 
investigations of full general relativity [^]. 

Finally, in the main text we have restricted ourselves to vacuum equations. The addition 
of a massless scalar field is straightforward because the Hamiltonian framework with which 
we began can be easily extended to accommodate a scalar field [3] • The main features that 
are generally expected from the analysis of Andersson and Rendall [3| are reflected in the 
resulting truncated system. Thus, if the energy density in the scalar field is small, one again 
has Bianchi II transitions and spikes. However, once the energy density exceeds a critical 
value, these disappear and the asymptotic dynamics at any spatial point is described just 
by the Bianchi model with a scalar field without any transitions. 

Derivations of analytical results and numerical simulations will appear in detailed papers 
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